clear;
  load('datamarketG2')


dt=0.01;
T1=20;
t1=0:dt:T1;
nt1=length(t1);
nmk(1)=ss(1);
ds=ss(2)-ss(1);


gnmk2=phiXmk2(1:flagstarmk2)-phiImk2(1:flagstarmk2);


for i=2:nt1
nmk(i)=nmk(i-1)+gnmk2(floor((nmk(i-1)-ss(1))/ds)+1)*dt;
end
qx2=1./(1-etaI*ixx);
  for i=1:nt1
rxmkt2(i)=rx2(floor((nmk(i)-ss(1))/ds)+1);
rpxmkt2(i)=rpx2(floor((nmk(i)-ss(1))/ds)+1);
qxmkt2(i)=qx2(floor((nmk(i)-ss(1))/ds)+1);
  end
  
  
  
 subplot(1,3,1)

plot(t1,rxmkt2,'LineWidth',1.5)
 hold on;

xlabel('$t$','interpreter','latex','FontSize',12)
title('A. $r_{t}^f$','interpreter','latex','FontSize',12)

axis('square')
  subplot(1,3,2)

plot(t1,rpxmkt2,'LineWidth',1.5)
 hold on;

 
 axis('square')
 
xlabel('$t$','interpreter','latex','FontSize',12)
title('B. $rp_t$','interpreter','latex','FontSize',12)

   subplot(1,3,3)

plot(t1,qxmkt2,'LineWidth',1.5)
 hold on;

 axis('square')
 
xlabel('$t$','interpreter','latex','FontSize',12)
title('C. ${\bf q}_t$','interpreter','latex','FontSize',12)




print -depsc figure_base_rrpq_stateGt5_new
 
 